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Proinflammatory responses by Toll-like receptors (TLRs) to malaria infec- 
tion are considered to be a significant factor in suppressing pathogen 
growth and in disease control. The key protozoan parasite Plasmodium fal- 
ciparum causes malaria through glycosylphosphatidyUnositols (GPIs), 
which induce the host immune response mainly via TLR2 signalhng. 
Experimental studies have suggested that malarial GPIs from P. falciparum 
are recognized by the TLR2 subfamily. However, the interaction site and 
their involvement in the activation mechanism are still unknown. A better 
understanding of the detailed structure of the TLR-GPI interaction is 
important for the design of more effective anti-malarial therapeutics. We 
used a molecular docking method to predict the binding regions of malarial 
GPIs with the TLR2 subfamily members. We also employed molecular 
dynamics simulations and principal component analysis to understand 
Kgand-induced conformational changes of the TLR2 subfamily. We 
observed the expected structural changes upon ligand binding, and signifi- 
cant movements were found in loop regions located in the ligand-binding 
site of the TLR2 subfamily. We further propose that the binding modes of 
malarial GPIs are similar to lipopeptides, and that the lipid portions of the 
Hgands could play an essential role in selective dimerization of the TLR2 
subfamily. 



Introduction 

Malaria is a predominant disease and is one of the 
most deadly infectious diseases [1]. In 2010, an esti- 
mated 3.3 bilhon people were in danger of contracting 
malaria and 655 000 deaths as a result of malaria were 
reported in 106 malaria-endemic countries [2]. The 
main cause of severe malaria include renal failure, 
pulmonary dysfunction, shock, disseminated intravas- 
cular coagulation, hypoglycaemia, metabolic acidosis 
and fiver dysfunction [3]. In humans, malaria is caused 
mainly by five protozoan parasite species of the 
genus Plasmodium: Plasmodium falciparum, Plasmo- 



dium vivax, Plasmodium malariae, Plasmodium ovale 
and Plasmodium knowlesi [4,5]. Plasmodium falciparum 
is the most deadly among these five species and con- 
tributes to 80-90% of the total malarial mortafity rate 
worldwide [6]. 

Innate immune mechanisms play a significant role in 
preventing diseases, including malaria, and also initiate 
long-lasting acquired immunity to eliminate pathogens 
and generate immunological memory in the late phase 
of infection [5]. Glycosylphosphatidyhnositols (GPIs) 
of P. falciparum are a crucial component among 
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malarial toxins, and they induce the production of 
proinflammatory cytokines, including tumour-necrosis 
factor-a, interleukin (IL)-l, IL-12 and nitric oxide, 
which are implicated in malaria pathogenesis [7]. Toll- 
like receptors (TLRs) are the most studied among the 
innate immune receptors, and it has been reported that 
TLRs, most notably TLR2, are involved in counteract- 
ing P. falciparum by recognition of GPIs, the anchor- 
ing protein of the parasite [8,9]. 

GPIs are glycoUpids that are pervasive in eukaryotes 
such as protozoa, fungi, yeast, plants, mollusks, 
insects and vertebrates [10]. The structure and function 
of GPIs vary from species to species, although the 
main function of GPIs in parasites is to anchor pro- 
teins to the cell membrane [6]. For P. falciparum, 
intact GPI or Man4-GPI, have four mannose (Man) 
groups, one glucosamine and one inositol (Ino); the 
tetramannosylglucosaminyl pentasaccharide is linked 
to phosphatidylinositol (Fig. SI) [6]. The Ino residue 
is acylated at the c-2 position. There are three fatty 
acid substituents in the Ino residue; each of the sn-1 
and c-2 positions have a saturated fatty-acyl substitu- 
ent and the sn-2 position has an unsaturated fatty-acyl 
moiety [6,10-12], 

TLRs are type 1 transmembrane proteins with an 
extracellular domain (ECD) for ligand recognition at 
the A'-terminus and a transmembrane helix domain, 
which connects the cytoplasmic Toll/IL-IR homology 
(TIR) domain that initiates signalhng cascades [13]. A 
molecular evolution study clustered vertebrate TLRs 
into six subfamilies: TLRl/2/6/10, TLRS, TLR4, 
TLRS, TLR7/8/9 and TLRl 1/12/13/21/22/23 [14]. The 
ECDs of these six subfamiHes of TLRs each bind to 
specific pathogen-associated molecular patterns. There 
are also possibilities of recognition of ligands by TLRs 
specific to species [15], Each ECD comprises 18-25 
copies of leucine rich repeats (LRRs) that form a 
horseshoe-like structure and are involved in various 
physiological functions, including immune responses 
and signal transduction [16]. To date, the ECD dimeric 
structures of TLRs 1, 2 and 4 (human), TLRs 2, 3 and 
6 (mouse), and TLRS (zebrafish) with their hgands 
have been reported [17-21]. These pubhshed TLR- 
ligand structures clearly exhibit the mechanism by 
which pathogen-associated molecular patterns are rec- 
ognized by their receptors [22]. In addition, these 
reported structures indicate that TLRs 1, 2, 4 and 6 
arc members of the 'atypical' subfamily of the LRR 
superfamily, whereas TLR3 and S are 'typical' subfam- 
ily members of the LRR superfamily [16,18]. Struc- 
tural differences between these two subfamilies in their 
LRR domains arise from the irregular arrangements of 
|3-sheets; the 'atypical' LRR subfamily can be divided 
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into A^-terminal, central and C-terminal subdomains 
[16]. The central domain is predicted to play a key role 
in ligand recognition [16]. 

The ECD of TLR2 can sense ligands from various 
pathogens, including parasites such as bacteria, fungus 
and viruses, and ligands such as triacylated peptides, 
lipoteichoic acid, diacylated lipoproteins, GPIs, phosp- 
holipomannan and glycoprotein B [23-27]. Homology 
models of TLRIO complexes were built and refined 
through molecular dynamics (MD) simulations to pre- 
dict the protein-ligand complexes hTLR2-hTLR10- 
Pam3CSK4, hTLRl-hTLR10-PamCysPamSK4 and 
hTLR10-hTLR10-PamCysPamSK4 [28]. X-ray crystal- 
lography studies of dimeric complexes of TLR2 in 
combination with TLRl, which recognizes Pam3CSK4, 
and another combination of TLR2 with TLR6, which 
recognizes Pam2CSK4, have dissected ECDs of the 
TLR2 subfamily [17,18]. Both of these dimeric crystal- 
lography studies revealed the structural basis for 
Hgand-induced receptor dimerization and explained the 
requirement for TLR2-TLR1 and TLR2-TLR6 hetero- 
dimerization in response to triacylated and diacylated 
lipopeptides, respectively, to initiate TLR signalling. 

Previous biochemical studies are available for GPIs 
recognition by TLR2 that leads to the activation of 
downstream signalling. The recognition of GPIs of the 
protozoans P. falciparum or Toxoplasma gondii 
appears to be via TLR2 and TLR4 [29]. In an experi- 
mental study by Krishnegowda et al. [30], using mouse 
macrophages and human monocytes, P. falciparum 
malarial GPIs consisting of three fatty acid chains 
were favourably recognized by human and mouse 
TLR2-TLR1 [30]. Moreover, one of the derivatives of 
GPIs called sn-2-\ys,o GPI was the ligand for the 
hTLR2-hTLR6 complex. The above result was con- 
firmed in another recent experimental study using mac- 
rophages from gene knockout mice, in addition to 
human monocytes and anti-human TLRl and TLR6 
sera [31]. The ECD of TLR2 has the potential to rec- 
ognize GPIs in the same binding sites of Upopeptides 
because the structural patterns of GPIs and lipopro- 
teins are similar, although they are different classes of 
compounds [30]. There is sufficient evidence for TLR2 
recognition of GPIs; however, the binding site of GPIs 
and the interacting residues in the protein that would 
be useful for developing anti-malarial drugs or vac- 
cines are still unknown. 

In the present study, we used some of the methods 
discussed below to determine the details of the 
interaction of the TLR2 subfamily with P. falciparum 
Man4-GPI and the sn-2 lyso GPI derivative. Molecular 
docking is a widely used modeUing tool for predicting 
the exact positioning of a ligand in the active site of a 
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protein [32]. Hence, in the present study, we employed 
molecular docking to investigate the interactions 
between P. falciparum Man4-GPI and hTLR2-hTLRl 
and between sn-2 lyso GPI and mTLR2-mTLR6. In 
addition, MD simulations that can report at the 
atomic level are appropriate for highlighting the 
dynamics of a given structure to validate the experi- 
mental studies on the ligand-induced dimerization 
analysis of TLRs [33]. It is well known that ligands 
induce dimerization of the TLR2 subfamily [17,18]; 
therefore, by utilizing MD techniques, we simulated 
the subfamily of TLR2 for 15 ns as a monomer and 
dimer in the absence and presence of the GPI to better 
understand the ligand-induced dimerization and activa- 
tion mechanism at the atomic level. 

Finally, we were interested in determining the most 
dominant or functionally important movements of 
the proteins when the dimers were bound and 
unbound to the hgands. Principal component analysis 
(PCA) was used to identify geometrical differences 
and global motions of the TLR2 subfamily when 
bound and unbound to GPIs. We observed that the 
Man4-GPI with three fatty acid moieties was recog- 
nized by hTLR2-hTLRl, and sn-2 lyso GPI with two 
fatty acid moieties was recognized by mTLR2- 
mTLR6. Hence, our observations are consistent with 
the hypothesis that the lipid moiety of the ligand 
plays a crucial role in TLR2 subfamily specificity of 
hgand recognition [31]. To our knowledge, no prior 
crystallography studies or molecular modelling studies 
have reported the interaction of GPIs with TLRs. 
The present study is the first to present the structural 
interactions of the GPIs with TLR2 subfamily mem- 
bers; these data could be used to facilitate the devel- 
opment of anti-malarial agents that specifically target 
TLR2 signalling. 



Results and Discussion 

Interaction between GPIs and TLR2 
heterodimeric complexes 

Although every TLR-ligand interaction is different, 
they all produce an 'm'-shaped dimeric complex with 
the C-terminus in the centre and the iV-terminus on 
the outside [34]. Among the TLRs, TLR2 has a unique 
characteristic in that it can form heterodimers with 
TLRl or TLR6 to recognize ligands [8]. Many experi- 
mental studies have suggested that Man4-GPI and sn-2 
lyso GPI are recognized by TLR2-TLR1 and TLR2- 
TLR6, respectively, and their interactions could be 
analogous to the recognition of lipopeptides [30,31]. 
From previous experimental studies [30,31], it is well 
known that GPIs bind TLR2 heterodimers; therefore, 
we predicted the binding sites of Man4-GPI and sn-2 
lyso GPI in the hTLR2-hTLRl and mTLR2-mTLR6 
complexes, respectively. Docking studies were per- 
formed using GOLD and mvd software to obtain insight 
into the binding conformation of Man4-GPI and the 
sn-2 lyso GPI derivative with TLR2 subfamily mem- 
bers. A total of 100 different conformations were 
obtained for both complexes; the top ten poses were 
selected for each complex, and their conformational 
stability was assessed. Based on the binding mode 
analysis, the optimal conformation for both Man4-GPI 
and sn-2 lyso GPI was selected from GOLD yield 
complex. In addition, the selected binding modes were 
similar to the lipopeptides conformation reported pre- 
viously [30,31]. 

Analysis of the selected hTLR2-hTLRl-Man4-GPI 
complex revealed that Man4-GPI was recognized in 
the convex region, which is in the same site where 
Pani3CSK4 binds [17] (Fig. 1). Man4-GPI has three 
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Fig. 1. Docked structure of the hTLR2-hTLR1-Man4-GPI complex. Ribbon diagram of hTLR2-hTLR1-Man4-GPI (hTLR2, hTLRI and l\/lan4-GPI 
are shown in green, blue and dark red, respectively). (A) Side view. (B) Top view. LRR modules and helix numbers are written in black. 
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lipid chains (Fig. SI) and, among those, the acyl chain 
at the sn-1 position and Ino predominate chain were 
buried in the hTLR2 pocket (volume: approximately 
1600 A^). The fatty acid substituent at the sn-2 posi- 
tion was accommodated into the TLRl channel (vol- 
ume: approximately 400 A^), and the head group of 
Man4-GPI was exposed out of the TLR2 pocket and 
interacted with the hydrophiUc residues (Fig. 2A). Six 
hydrogen-bonding contacts between hTLR2-hTLRl 
and the Maru-GPI were observed (Fig. 2B), and the 
TLR residues that participated were Phe325, Asp327, 
Asn294 of hTLR2, as well as Gly313 and Gly316 of 
hTLRl. 

The best docking complex for mTLR2-mTLR6-5«-2 
lyso GPI among the top ten ranked complexes was 
selected for analysis. As stated above, the orientation 
resembles the available crystal structure of the 
mTLR2-mTLR6-Pam2CSK4 complex [18]. The volume 
of mTLR2 pocket is approximately 1200 A^, although 
the hydrophobic channel of mTLR6 was blocked by 
phenylalanine residues 343 and 365 [18]; however, 
mTLR2-mTLR6 could still recognize sn-2 lyso GPI 
(Fig. 3) because it lacks the fatty acid substituent at 
the sn-2 position of Man4-GPI (Fig. SI). As a result of 
the presence of the two phenylalanine residues present 
in mTLR6, we failed to produce Man4-GPI in a com- 
plex with mTLR2-mTLR6. Eight hydrogen bonds 
were observed in the mTLR2-mTLR6-OT-2 lyso GPI 
complex (Fig. 4). The amino acids that formed hydro- 
gen bonds with sn-2 lyso GPI were Glu264, Asp294, 
Asp327 of mTLR2, and Phe319 and Lys321 of 
mTLR6. Asp327 was the common hydrogen-bonding 
residue in the docked complexes (hTLR2-hTLRl- 
Man4-GPI and mTLR2-mTLR6-i-«-2 lyso GPI), as 
observed in the previously reported hTLR2-hTLRl- 
Pam3CSK4 and mTLR2-mTLR6-Pam2CSK4 crystal 
structures. The hTLR2 contains an Asn residue, 
whereas mTLR2 contains an Asp residue at position 
294, and both of these residues were involved in 
hydrogen bonding with their respective hgands. In 
addition to crystallographic studies on TLR2 [17,18], 
there is also a recent mutational study that supports 
the docked binding site [35]. Phe325 formed a hydro- 
gen bond with Maru-GPI, and Asp327 formed a 
hydrogen bond with Man4-GPI and sn-2 lyso GPI; 
these residues were shown to be necessary for mediat- 
ing cellular responses and for strong nuclear factor-KB 
activation in TLR signalling [35]. In another mutagen- 
esis study conducted on TLRl, the Gln316Lys substi- 
tution was shown to abrogate hgand stimulation, 
indicating its importance for Hgand binding [36]. The 
docked hTLR2-hTLRl-Man4-GPI complex showed a 
potential hydrogen bond between Gln316 of TLRl 



and the GPI, which is consistent with previously pub- 
hshed data and supports the validation of the selected 
docked complex. 

Structural comparison of GPI and lipopeptide 
complexes of the TLR2 subfamily 

TLRs usually recognize ligands at their concave sur- 
face, although crystallographic studies have revealed 
that, when TLR2 heterodimerizes with TLRl or 
TLR6, ligand recognition occurs at their convex region 
[34]. The selected docked complexes of Man4-GPI and 
sn-2 lyso GPI also confirmed this unique characteristic 
of TLR2 (Figs 2 A and 4A). Moreover, the convex 
region ligand-binding pocket of TLR2 that heterodi- 
merizes with TLRl and TLR6 to recognize the 
Hpopeptides is present at the edge of the central and 
C-terminal domains [16]. The Man4-GPI and sn-2 lyso 
GPI binding sites were the same as observed for 
Pam3CSK4 and Pam2CSK4, respectively. Interestingly, 
some of the interacting residues in the crystal structure 
of the hTLR2-hTLRl-Pam3CSK4 complex were also 
observed to participate in the docked hTLR2-hTLRl- 
Man4-GPI complex (Fig. 5 and Table 1). To obtain a 
better view of the conserved residues in both com- 
plexes, we have compared them in Fig. S2. In particu- 
lar, the hTLR2 residues (Asn294 and Asp327) and 
hTLRl residues (Gly313 and Gly316) participate in 
the hydrogen-bonding network within the Pam3CSK4 
crystal complex and the predicted Man4-GPI complex. 
In the crystal structure, two ester-bound hpid chains 
and a single amide-bound lipid chain from Pam3CSK4 
bound to the hTLR2 and hTLRl pockets, respectively. 
Similarly, sn-l and c-2 fatty acid chains shared the 
hTLR2 pocket, and the sn-2 fatty acid chain bound to 
hTLRl. The glycerol and peptide portions in the head 
group of Pam3CSK4 formed hydrogen bonds and 
hydrophobic interactions wdth hTLR2-hTLRl. Simi- 
larly, analysis of the docked hTLR2-hTLRl-Man4- 
GPI complex indicated hydrogen-bonding interactions 
between the head group residues of Man4-GPI (Man, 
glucosamine and Ino) and the TLRs (Fig. 2). 

A notable structural similarity of sn-2 lyso GPI and 
Pam2CSK4 is the number of lipid chains, which allows 
sn-2 lyso GPI to use the Pam2CSK4 binding pocket in 
the mTLR2-mTLR6 heterodimer. Therefore, both 
hgands have few commonly interacting residues (Fig. 5 
and Table 1; see also Fig. S3). Pam2CSK4 has two 
ester-bound lipid chains and sn-2 lyso GPI also has 
two fatty acid chains at the sn-1 and C-2 position of 
the Ino moiety. Both available fatty acid moieties in 
the sn-1 and c-2 positions of Ino of sn-2 lyso GPI were 
bound in the mTLR2 pocket, which is the same region 
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Man.-GPI 




Fig. 2. The Man4-GPI binding site in tlie docked structure of the lnTLR2-lnTLR1 complex. (A) Tine lnTLR2 and InTLRI residues involved in 
Man4-GPI binding are shown in blue and green, respectively. The hydrogen bonds are shown by dotted red lines. Carbons, nitrogens, 
oxygens and phosphorous of Man4-GPI are coloured in orange, blue, red and purple, respectively. (B) The chemical structure of Man4-GPI 
and residues forming hydrogen bonds. Hydrogen bonds are shown with dotted red lines, and their bond distances are written above the 
lines in angstroms. Carbons, nitrogens, oxygens and phosphorous of Man4-GPI are coloured in black, blue, red and purple, respectively. 
Apostrophes are used for hTLRI residues to differentiate from residues in hTLR2. 
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Fig. 3. Docked structure of the mTLR2-mTLR6-sn-2 lyso GPI complex. Ribbon diagram of mTLR2-mTLR6-sn-2 lyso GPI (mTLR2 and mTLR6 
witli sn-2 lyso GPI are slnown in cyan and grey, respectively). (A) Side view. (B) Top view. LRR modules and helix numbers are written in black. 



where Pam2CSK4 lipid chains are bound [18]. There- 
fore, the small size of the mTLR6 pocket does not 
make a significant impact on .m-2 lyso GPI binding. 
The Phe319 in the LRRU loop of mTLR6 forms 
hydrogen bonds with .m-2 lyso GPI and Pam2CSK4; 
however, Phe314, the corresponding residue in hTLRl, 
does not appear to interact with the ligands, which 
may be a result of differences in the heterodimerization 
of TLRs. 

To evaluate the docked complexes of hTLR2- 
hTLRl-Man4-GPI and mTLR2-mTLR6-.v«2 lyso 
GPI, we predicted the binding affinity using x-score 
[37] and compared it with the crystal structures of tri- 
and di-acyl lipopeptides. Table 2 shows the calculated 
binding energy of hTLR2-hTLRl-Man4-GPI value to 
be —13.4 kcal-moP', which is close to the value of 
— 12 kcal moP' obtained for the crystal structure of 
hTLR2-hTLRl-Pam3CSK4. Similarly, the mTLR2- 
mTLR6-.y«-2 lyso GPI had a calculated binding energy 
of —12 kcal-moP', which is almost same as the bind- 
ing energy of —12.2 kcal-moP' predicted for the 
mTLR2-mTLR6-Pam2CSK4 complex. Slight differ- 
ences in the binding energy may be a result of differ- 
ences between the GPI and lipopeptide complexes such 
as: (a) P. falciparum GPIs have long lipid chains com- 
pared to the lipopeptides that strongly bind their 
respective receptors and (b) GPIs form more hydrogen 
bonds with TLRs than with lipopeptides. It should be 
noted that the crystal structures of TLR2 with other 
ligands such as pnLTA and PE-DTPA have shown a 
lack of hydrogen bonds with TLR2, which leads to 
low binding affinity because they are unable to stimu- 
late dimerization or activate the receptor [18]. 

TLR2 activates the same signalling pathway irre- 
spective of whether it heterodimerizes with TLRl or 
TLR6 [38]. Because our results show that malarial 



GPIs bind with TLR2 in complex with TLRl or 
TLR6, they resemble the known hpopeptides com- 
plexes [17,18]. Nevertheless, the P. falciparum GPIs 
and lipopeptides have a similar series of lipid chains; 
hence, these P. falciparum GPIs can possibly induce 
similar TLR2 signalling induced by hpopeptides. 
Although it has been reported that GPIs from differ- 
ent species with different structures are diversely recog- 
nized by TLRs, similar to lipopeptides, after the GPIs 
bind to LCDs of TLR2, they may induce the dimeriza- 
tion of LCDs that results in dimerization of intracellu- 
lar TIR domains, which leads to the recruitment of 
adapter molecules to initiate signalhng [16]. The GPI- 
induced TLR2 signalling, detailed with the adapter 
molecules, is shown in Fig. 6. 

Structural stability analysis from MD simulations 

Protein stability during MD simulations was moni- 
tored by RMSD to examine differences in the mono- 
meric and dimeric forms of the TLR2 subfamily, as 
well as their effect in the presence of ligand. Therefore, 
we included dimeric complexes without hgands, even 
though it is not experimentally possible in the TLR2 
subfamily. Accordingly, 15-ns MD simulations were 
carried out for two sets of TLRs. The first set con- 
tained four systems: hTLR2, hTLRl, hTLR2-hTLRl 
and hTLR2-hTLRl-Man4-GPI; the second set con- 
sisted of four systems: mTLR2, mTLR6, mTLR2- 
mTLR6 and mTLR2-mTLR6-OT-2 lyso GPI. In addi- 
tion, to provide a control to compare RMSD of the 
docked complexes, the 15-ns MD simulations for the 
crystal structure of hTLR2-hTLRl-Pam3CSK4 were 
carried out (Fig. S4). The plot of the RMSD of all 
backbone atoms as a function of time revealed relative 
changes in the structure, indicating significant protein 
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Fig. 4. The sn-2 lyso GPI binding site in the docked structure of the mTLR2-mTLR6 complex. (A) The mTLR2 and mTLR6 residues involved 
in sn-2 lyso GPI binding are shown in cyan and grey, respectively. The hydrogen bonds are shown by dotted red lines. Carbons, nitrogens, 
oxygens and phosphorous of the sn-2 lyso GPI are coloured in yellow, blue, red and purple, respectively. (B) The chemical structure of sn-2 
lyso GPI and residues forming hydrogen bonds. Hydrogen bonds are shown with dotted red lines, and their bond distances are written 
above the lines in angstroms. Carbons, nitrogens, oxygens and phosphorous of sn-2 lyso GPI are coloured in black, blue, red and purple, 
respectively. Apostrophes are used for mTLR6 residues to differentiate from residues in mTLR2. Residue Asp327 forms multiple hydrogen 
bonds with the ligand and so it is shown more than once. 



mobility. As shown in Fig. 7A, the first set of RMSD 
values for monomers were within a fluctuation range 
of 0.9-3.3 A throughout the simulations, except the 
dimers. For hTLR2-hTLRl-Man4-GPI, the RMSD 
values were in the range 0.9^.7 A, whereas, for the 



hTLR2-hTLRl system, the RMSD values were in the 
range 0.9-5.1 A during the 15-ns trajectory. Accord- 
ingly, within 15 ns, the RMSD of the hgand-bound 
hTLR2-hTLRl-Man4-GPI system tended to converge, 
and a longer simulation time may be needed to reach 
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Central domain 

LRR9 xLxxLxLx B2 

niTLR2 VSSPMKKLAFRGSvItDEsIneBklIR 274 

hTLR2 TNSLIKKFTFRNVKITDESLFQvJkLLN 274 

*AA/\AAAAAA/\ H2' 

mTLR6 RGPTLLNVTLQHIETTWKCSVKLFQFF 273 

hTLRl TNPKLSSLTLNNIETTWNSFIRILQLV 268 

LRRIO xLxxLxLx 

[nTLR2 YILELSE|e|DDCT|NGLGD|NPSESDWS 304 

hTLR2 QISGLLELEFDDCTLNGVgIfRASDNDRVI 304 

mTLR6 WPRPVEYLNIYNLTITERIDREEFTY 299 

hTLRl WHTTVWYSSISNVKLQGQLDFRDFDY 294 

UtRll xLxxLxLx 

mTLR2 e|GKVET|t|RR|h|PQFYL|yBst|YS 333 

hTLR2 DPGKVETLTIRR§HIPRFYliY|Di,STLYS 333 

*A/\A/W\ H3' 

mTLR6 SETALKSLMIEHVKNQVfBsKEALYSVF 328 

hTLRl SGTSLKAI.s|H0WSDVF||FlQ5YiyEl| 323 

C- terminal domain 

LRR12 xLxxLxLxxN- w^wv/v H3 

mTLR2 Bs'<v:-;ftITVF.i:.-:K|FLV|is|sc 357 

hTLR2 |ter|kbJ|tvei;skvfl|pclJsq 357 

mTLR6 AEMNIKMLSISDTPFIHMVCPP 350 

hTLRl SNMNIKNFTVSGTRfVHMLCPS 3<5 



I Mai\4-GPI and Pam,CSK4 binding iPam^CSK* and sn-2 lyso GPI binding 
[] MaOi-GPI binding sn-2 lyso GPI binding 

I PamjCSK. b«K)ing | PamjCSK. binding 

Fig. 5. Structure-based sequence alignment of the ligand-binding 
regions of hTLR2, hTLRl, mTLR2 and mTLR6. The ligand binding 
regions of hTLR2, hTLRl, mTLR2 and mTLR6 sequences are 
aligned based on their structures. Conserved leucines and residues 
in the asparagine ladder are written above the sequences. The 
positions of a-helices are indicated by coils above the sequences 
and are numbered. 

relaxed equilibrium state. Monomers and the ligand- 
unbound liTLR2-liTLRl system had not reached equi- 
hbrium up to 15 ns. The RMSD determined for the 
hgand-unbound hTLR2-hTLRl was larger (5.2 A) 
than the value calculated for the ligand-bound 
hTLR2-hTLRl-Man4-GPI complex (4.7 A). These 
results were expected because it has been previously 
shown that both monomers and hgand-unbound 
dimers are known to be inactive [39]. These data sug- 
gest that hgands can induce/stabilize substantial con- 
formation changes on the dimeric complex, whereas 
the ligand-free dimer was less stable during the MD 
simulations. Analysis of the second set showed that 
the structural drift was continual in all the structures, 
and the RMSD was in the range 0.9-5 A, except for the 
ligand-bound niTLR2-mTLR6-.y«-2 lyso GPI complex 
(Fig. 7B). The RMSD value of the mTLR2-mTLR6- 
sn-2 lyso GPI was in the range 0.9-3.2 A, although, 



Table 1. Hydrogen-bonding interactions of members of the TLR2 
subfamily with malarial GPIs and lipopeptides. The hydrogen- 
bonding interactions of l\/lan4-GPI, Pam3CSK4, sn-2 lyso GPI and 
Pam2CSK4 with their respective TLRs and hydrogen bond 
distances are shown. 





Residue name 




Hydrogen bond 


TLRs 


and number 


Ligand name 


distance (A) 


hTLR2 


Asn294 


Man4-GPI 


2.7 






Pam3CSK4 


3 




Phe325 


Man4-GPI 


2.9 




Asp327 


Man4-GPI 


2.6 and 2.5 






Pam3CSK4 


2.9 




Phe349 


Pam3CSK4 


3 


hTLRl 


GlySia 


Man4-GPI 


2.6 






Pam3CSK4 


2.9 




Gln316 


Man4-GPI 


2.8 






Pam3CSK4 


2.6 


mTLR2 


Glu264 


sn-2 lyso GPI 


2.7 




Asp294 


sn-2 lyso GPI 


2.7 




Phe325 


Pam2CSK4 


3.2 




Asp327 


sn-2 lyso GPI 


2.5, 2.5, 2.5 








and 2.6 






Pam2CSK4 


2.5 




Phe349 


Pam2CSK4 


3 


mTLR6 


Phe319 


sn-2 lyso GPI 


2.9 






Pam2CSK4 


2.9 




Lys321 


sn-2 lyso GPI 


2.6 



after 3 ns, the system was equilibrated and became sta- 
ble throughout the simulations and reached a final 
value 3.1 A. In both sets, the results showed lower 
mean RMSD values for hgand-bound forms than in 
the Kgand-free forms (3.3 A for hTLR2-hTLRl- 
Man4-GPI, 3.7 A for hTLR2-hTLRl, 2.7 A for 
mTLR2-mTLR6-«-2 lyso GPI, and 2.9 A for mTLR2- 
mTLR6; Fig. 7), indicating that binding of ligand to 
this receptor induced/stabilized complex formation. An 
RMSD analysis of the crystal structure of hTLR2- 
hTLRl-Pam3CSK4 indicates that the RMSD was in 
the range 0.9^ A and it is more stable compared to 
experimentally solved complexes of hTLR2-hTLRl- 
Man4-GPI. Hence, the results show that the hTLR2- 
hTLRl-Man4-GPI complex needs more time to 
stabilize compared to the control. It should be noted 
that, in a MD simulation study on myeloid differentia- 
tion protein-lipid A-TLR4 complex, the presence of 
ligand energetically stabilized the complex [39]. To 
observe the conformational changes on the TLR2 sub- 
family during MD simulation, the distances between 
the A'-terminal and C-terminal ends of the TLR2 sub- 
family structures were calculated for the initial struc- 
tures and after MD runs using pymol, version 1.5.0.4 
(Schrodinger, LLC, Portland, OR, USA; Table 3). 
Interestingly, after MD simulations, the distances 
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Table 2. Binding affinity prediction. Tiie mean binding affinity scores of docked protein-ligand complexes of hTLR2-hTLR1-l\/1an4-GPI and 
mTLR2-mTLR6-sn-2 lyso GPI and crystal structures of liTLR2-hTLR1-Pam3CSK4 and mTLR2-mTLR6-Pam2CSK4. pKa, negative logarithm of 
dissociation constants. 





Hydrophobic pair 


Hydrophobic match 


Hydrophobic surface 


Predicted mean 


Predicted binding 


Protein-ligand complex 


score ipKa) 


score (p^Cd) 


score ipKa) 


binding affinity ipK^) 


energy (kcal-mor^) 


lnTLR2-hTLR1-Man4-GPI 


10.91 


10.31 


8.48 


9.9 


-13.4 


hTLR2-hTLR1-Pam3CSK4 


9.24 


10.23 


7.12 


8.86 


-12 


mTLR2-mTLR6-sn-2 lyso GPI 


10.03 


9.21 


7.39 


8.88 


-12 


mTLR2-mTLR6-Pam2CSK4 


9.97 


9.83 


7.18 


9 


-12.2 



'^y hTLR2-hTLR1-»««n.-GPl 



mTlJ?2-mTLR6-«/^ 




K(K comptox 
(IKK-v. IKK-oft 
IKK-«) 



0«Mh domain of MyDBS. 

IRAK4-1'2 complex 



MARK (ERK. JNK. p38)-TAK1-TAB 
TRAF6-UBC13*UEV1A 




TAK1 




Ap-1 



hl-B-NF-«B (pSO-peS) 





Induce the expression o( ^ 

♦*** pfoinnemmalory genes ^ w * ' ' ' '.^"""^"'^ 

NFHcBMndtloONA Ap-I btnd. to 0«A 

»•*'* Nuclear r 

'" 



Fig. 6. Schematic view of 
Plasmodium falciparum GPI-induced TLR2 
signalling. TLR2 ECD recognizes GPIs 
either in combination with TLRl or TLR6, 
then their TIR domains dimerize on the 
cytosolic face and undergo structural 
changes to recruit adaptor molecules: 
MyD88, MAL, IL-1 R-associated kinases 
(IRAKs), (TGF-p)-activated kinase (TAK1), 
TAKI-binding protein 1 (TABl) and TAB2, 
and then tumour-necrosis factor-receptor- 
associated factor 6 (TRAF-6) forms a 
complex with ligases Ubc13 and UevIA, 
which leads to the activation of the 
mitogen-activated protein kinase and 
nuclear factor-KB pathways to induce the 
production of proinflammatory cytokines. 



between both their terminals were decreased from their 
initial structures when hTLR2 and mTLR2 were 
monomers or dimers with hgands. TLR2 is the key 
receptor in the heterodimers and the reduced distance 
observed between both the ends is likely to enable 



ligand binding and dimerization with TLRl or TLR6. 
By contrast, for hTLRl and mTLR6, the distance 
between the two ends was found to be increased, 
which may be to ease ligand binding and heterodimer- 
ize with TLR2. The distances measured between the 
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A 0.6 



0.5 



0.4 



RMSD 



hTLR2 
hTLR1 



hTLR2-hTLR1 
hTLR2-hTLR1-Man,-GPI 




Time (ns) 



B 0.6 



RMSD 




Time (ns) 



Fig. 7. RMSD plots of ligand-bound and ligand-unbound complexes. (A) The RMSD plot for the backbone atoms of hTLR2, hTLRI, hTLR2- 
hlLRI and hTLR2-hTLR1-Man4-GPI relative to the initial structure over 15-ns MD simulations. (B) The RMSD plot for the backbone atoms of 
mTLR2, mTLR6, mTLR2-mTLR6 and mTLR2-mTLR6-sn-2 lyso GPI relative to the initial structure over 15-ns MD simulations. 



Table 3. Comparison of TLR2 subfamily ECDs before and after MD simulations. The distances between the N- and C-terminal ends of TLR2 
subfamily members were calculated for their initial and 15-ns structures. 



TLRs 


Initial structure 
(A) 


After MD simulations 
as monomer (A) 


After MD simulations 

as ligand-unbound dimer (A) 


After MD simulations as 
ligand-bound dimer (A) 


hTLR2 


59.1 


49.3 


50.9 


52 


hTLRI 


44.1 


55.7 


43.3 


45 


mTLR2 


56.3 


52 


62.7 


53.9 


mTLR6 


47.9 


41.8 


48.4 


48.7 



A^-terminal and C-terminal ends for the final control 
MD simulation structures for hTLR2-liTLRl- 
Pam3CSK4 and the values for hTLR2 and liTLRl, are 
54.1 and 44.6 A, respectively. Those values are closer 
to the 52 A of hTLR2 and 45 A of hTLRI in the 
hTLR2-hTLRl-Man4-GPI complex (Table 3). In addi- 
tion, we have calculated the significant movements 
observed in the ligand-binding region (Table SI). In 
the case of hgand-unbound dimers, the values were var- 
iable and inconclusive because these complexes cannot 
be determined experimentally. Moreover, we calculated 
the distances of the hTLR2 and hTLRI both in ligand- 
bound and Ugand-unbound forms after 15-ns MD sim- 
ulations (Fig. S5). The distances for liganded and unli- 
ganded hTLR2 and hTLRI were 44.9 and 60.4 A, 
respectively. Interestingly, this difference shows that 
unliganded TLR2 heterodiniers were unstable and 
found to disassociate during the course of longer simu- 
lation. The values agreed with previously published val- 
ues [40,41], indicating that, upon Ugand binding, TLRs 
switch to a shortened and activated form. 



Residual flexibility analysis from MD simulations 

The final 5 ns were used to calculate the RMSF plot 
of backbone atoms for hTLR2, hTLRI, hTLR2- 
hTLRl, hTLR2-liTLRl-Man4-GPI, mTLR2, mTLR6, 
mTLR2-mTLR6 and mTLR2-mTLR6-.v«-2 lyso GPI 
(Fig. 8). In addition, the RMSF values were calculated 
for the crystal structure of hTLR2-hTLRl-Pam3CSK4 
to compare the fluctuations of residues in docked com- 
plexes (Fig. S6). For both ligand-unbound forms, the 
two sets showed large fluctuations for the backbone 
atoms compared to the corresponding ligand-bound 
form. Notably, the RMSF for control hTLR2- 
hTLRl-Pam3CSK4 also showed results similar to 
other docked complexes compared to ligand-unbound 
forms (Fig. S6). The results from both the crystal and 
docked complexes showed that the residues in the 
ligand-binding regions stabilized upon ligand recogni- 
tion. The RMSF values also indicated that the residues 
with higher fluctuation values were those located in 
the loop or helix regions. In both of the sets subjected 
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Fig. 8. RMSF Plots of ligand-bound and ligand-unbound complexes. RMSF plots of residue fluctuations obtained by averaging atomic 
fluctuations over the last 5 ns of MD simulations. Ligand-bound and ligand-unbound states of (A) hTLR2, (B) hTLRI, (C) mTLR2 and (D) 
mTLR6. 



to MD simulations, the A'-terminal and C-terminal 
regions showed higher fluctuation than other LRR 
regions. The TLR2 subfamily binds with a wide-range 
of microbial compounds and proteins in addition to 
hpoproteins and lipopeptides [25]. Among them, pepti- 
doglycan of Staphylococcus aureus can bind to the 
iV-terniinus or C-terniinus of the TLR2 [42]; proline- 
proline-glutamic acid 18 of Mycobacterium tuberculosis 
has been shown to interact with the C-terminal region 
of TLR2 [43]; and both the TV- and C-terminus LRRs 
of TLRl and TLR6 play a minimal role in accommo- 
dating LT-IIb-Bs [44]. Hence, we hypothesize that the 
N- and C-terminal fluctuations observed in our analy- 
sis are reasonable because they are involved in diverse 
interactions with multiple ligands. Moreover, the 



ligand-induced dimerization takes place in the C-termi- 
nal region of TLR2, TLRl and TLR6 [17,18]. In addi- 
tion to these A^- and C-terminal regions, there were 
high fluctuations on hTLR2, hTLRI and mTLR2 in 
their ligand-binding LRR 9-12 regions. For mTLR6, 
fluctuations were observed in the ligand-binding site 
LRR 11-14. 

In hTLR2 and hTLRI, relative differences in the 
loop regions were observed compared to the helix and 
P-sheet regions in the absence and presence of Man4- 
GPI. The binding site residues showed different flexi- 
bilities, demonstrating that these residues are stringent 
after binding to Man4-GPI (Fig. 8A,B). Interestingly, 
in the hgand-biiiding region, residues Glu246, Thr247, 
Asn248, Asn300, Asp301, Arg302 and Val303 of 



6206 



FEBS Journal 280 (2013) 6196-6212 © 2013 The Authors. FEES Journal published by John Wiley & Sons Ltd on behalf of FEBS 



p. Durai ef al. 



Binding modes of malarial GPIs in TLR2 subfamily 



hTLR2, as well as Asp288 and Phe289 of hTLRl, 
showed more fluctuations than the others in 
the absence and presence of Man4-GPI. The overall 
fluctuations of residues in the control hTLR2- 
hTLRl-Pam3CSK4 complex are the same as the 
hTLR2-hTLRl-Man4-GPI complex. Moreover, when 
we compared the residual fluctuations with individual 
monomers and unliganded dimers, we observed fewer 
fluctuations (Fig. S6). This observation indicated con- 
formational changes in these residues that were stabi- 
lized by an interaction with the Man4-GPI Ugand. 
Notably, a previous mutagenesis study on TLR2 
implies that the above stabilized residues Glu246 and 
Asp301 had an intermediate response to Pam3CSK4 
regardless of there being no contact with TLRl or the 
ligand [35]. In another set of mTLR2 and mTLR6, the 
RMSF plot showed that the binding of sn-2 lyso GPI 
to mTLR2-mTLR6 resulted in a decreased flexibility 
of the residues and was relatively stable (Fig. 8C,D). 
The fluctuations of residues in hTLR2 and mTLR2 
were not the same, possibly as a result of structural 
differences: (a) the lipid binding pocket of mTLR2 is 
smaller than hTLR2 and (b) sequence similarity 
between hTLR2 and mTLR2 is only 68%, which leads 
to structural changes between those TLRs in the 
LRRIO and LRRll regions [17]. As discussed above, 
the role and fluctuations of Glu246 in the hTLR2- 
hTLRl-Man4-GPI complex and the same residue pres- 
ent in the loop region before the beginning of LRR9 
of mTLR2 also showed high fluctuations. In addition, 
LRRIO residues, including Glu299 and SerBOO, showed 
fluctuations, although the RMSF value was relatively 
low for residues in the protein-ligand complex. The 
ligand-binding LRR 9-12 of mTLR2 and LRR 11-14 
of mTLR6 was also flexible. The RMSF plot clearly 
shows the overall increased fluctuations of ligand-inde- 
pendent complexes and the RMSF values for control 
agree with the results for the P. falciparum GPI com- 
plexes (Fig. S6). 

PCA 

To better understand the conformational changes of 
the hgand-induced dimerization mechanism of TLR2 
with their partners, 15-ns MD trajectories were ana- 
lyzed by examining only the motions with the principal 
components. The contributions to the motion of the 
first ten collective modes for all TLRs are shown in 
Fig. S7. The application of PCA to the backbone 
atoms of hTLR2, hTLRl, hTLR2-hTLRl and 
hTLR2-hTLRl-Man4-GPI motions of the simulation 
indicated that over 23%, 15%, 67% and 66%), respec- 
tively, of these motions are accounted for by the first 



two eigenvectors (Fig. S7). On the other hand, the first 
two eigenvectors covered ll%o, 16%, 12% and 40%) 
for mTLR2, mTLR6, mTLR2-mTLR6 and niTLR2- 
mTLR6-i'«-2 lyso GPI, respectively. We also investi- 
gated whether conformational changes occurred if the 
ligand was removed from the heterodimeric complexes 
of mTLR2 and hTLR2; therefore, we manually 
removed their Ugands and performed the PCA for 
their backbone atoms. PCA motions of the first two 
eigenvectors for hTLR2-hTLRl and mTLR2-mTLR6 
dimers without hgand were over 61% and 12%, 
respectively (Fig. S7). The motions of backbone atoms 
for the first two eigenvectors and their amplitude are 
given in Fig. S8. For dimers bound and unbound to 
their ligands, the first two modes covered the crucial 
fluctuations with motion. Consequently, the projec- 
tions of the TLR2 subfamily in each MD trajectory 
onto the essential space (planes) defined by eigenvector 
1 and 2 (Fig. 9) allow the visualization of the confor- 
mational spaces sampled from MD calculations. The 
distribution of all monomers of hTLR2, hTLRl, 
mTLR2 and mTLR6 were relatively similar, which 
indicates that, when TLRs remain as a monomer, they 
stay in a similar conformational state; however, the 
dimer with hgand showed dissimilar eigenvectors. The 
size of each cluster in Fig. 9 indicates that hTLR2- 
hTLRl-Man4-GPI and mTLR2-mTLR6-.v«-2 lyso 
GPI undergo large conformational changes during 
simulation. Each trajectory of the MD simulations 
projected onto eigenvectors 1 and 2 had two well- 
defined clusters, indicating that the monomer and 
dimer had two distinct minima during the molecular 
dynamics trajectory. The removal of ligand from the 
dimeric complexes of hTLR2-hTLRl and mTLR2- 
mTLR6 led to an extensive sampling of phase space 
until the system reached different local minima on the 
energy landscape. 

The protein domain motions were analyzed to iden- 
tify the corresponding motion modes by the first two 
eigenvectors. The 0-15-ns representative structures of 
both MD sets for the first two eigenvectors are shown 
in Fig. S9. As expected, the monomer of the TLRs 
showed large movements in the C-terminal region, as 
found in our RMSF results (Fig. 8). The C-terminal 
regions are responsible for interacting with other mole- 
cules because there were large conformational motion 
differences for the monomer compared to the dimer 
(ligand-bound form). These results are consistent with 
our RMSF results in that the C-terminal regions of 
TLR2 dimers fluctuated when they are bound to 
ligand, whereas the TLR2 dimeric complexes without 
Hgand showed even higher fluctuations in their C-ter- 
minal regions (Fig. 8). Furthermore, there were 
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Fig. 9. PCA of the TLR2 subfamily. The cloud represents the 15-ns trajectories projected onto the first two eigenvectors. Every fifth frame 
of the respective trajectories has been used in both projections. The x-axis and /-axis show the projection of the structures of the backbone 
atoms in the IVID trajectories onto the essential space (planes) defined by eigenvectors 1 and 2. (A) The clouds coloured in yellow, green, 
blue and red display the monomer of hTLR2, hTLRI, hTLR2-hTLR1 and hTLR2-hTLR1-l\/lan4-GPI trajectories, respectively. (B) The clouds 
coloured in yellow, green, blue and red display the mTLR2, mTLR6, mTLR2-mTLR6 and mTLR2-mTLR6- sn-2 lyso GPI trajectories, 
respectively. 



notable motions observed in the loops in and around 
the ligand-binding regions. The projection of eigenvec- 
tors 1 and 2 displayed a pivotal movement of loop 
regions, which distorted and exposed the basic concave 
surface of the TLR2 subfamily. The hTLR2 has a 
large concave cavity compared to that of mTLR2 
because its loop regions have more predominant 
motions than mTLR2. As can be seen in Fig. S6, the 
difference in motions of the LRRIO loop region of all 
forms of hTLR2 (hTLR2 as monomer, Man4-GPI- 
unbound hTLR2 dimer with hTLRI, and Man4-GPI- 
bound hTLR2 dimer with hTLRI) suggest that, during 
the MD simulations, the monomer of the liTLR2 loop 
region comes forward to accommodate the hgand to 
induce dimerization. Alternatively, Man4-GPI-bound 
hTLR2 showed that the LRRll loop region was 
shghtly extended. More fluctuations could be seen in 
LRR 9-11 loop regions for hTLR2 and loops in the 
LRRll for liTLRl. In the case of the monomeric 
form of mTLR2, motion of Hgand-binding loop 
regions was observed; however, the ligand-bound 
dimer of mTLR2 revealed that fluctuations in the 
LRR 9-11 loop region were completely stabilized by 
their interaction with Kgand. In particular, mTLR6 
showed large motions in their ligand-binding loop 
regions located in and around the LRRll loop, which 
is known to interact with lipopeptides. The monomer 
of mTLR6 showed large motions in ligand-binding 



loop regions, whereas the dimer with ligand showed 
less motion. The dimer of Kgand-unbound mTLR6 
showed no different motions in their ligand-binding 
loop regions and looked similar to monomer. The 
motions in the loop regions of monomeric mTLR6 
indicate that they provide space for hgands to bind 
with mTLR2. Subsequently, when mTLR6 interacts 
with ligand through dimerization with niTLR2, the 
loop regions were relatively fixed. Overall, all of these 
results showed that the transition of monomer 
involved bending motions of loops in the ligand-bind- 
ing regions, particularly in the LRR 9-1 1 modules. 
The observed changes in PCA are consistent with the 
RMSF analysis (Fig. 8), indicating that the loop 
regions in the ligand-binding central domain of TLR2 
undergo high fluctuations in the absence of ligand 
molecules to induce the dimerization with their sub- 
family members. 

The aim of the present study was to determine the 
binding region of P. falciparum malarial GPIs in the 
TLR2 subfamily, as well as to understand the ligand- 
induced dimerization mechanism at the atomic level. 
In the present study, molecular docking, MD simula- 
tions and PCA were employed to investigate the inter- 
action between malarial GPIs and TLR2 subfamily 
members. Conformational changes and interactions at 
their interfaces were analyzed. PCA suggested that the 
activation loops located in the central region of the 
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TLR2 subfamily were more flexible and displayed sig- 
nificant changes upon ligand binding. After the bind- 
ing of malarial GPIs, the A'- and C-terminal domains 
of TLR2 displayed a constrained conformation and 
the wiggling motion was localized to the activation 
loops. Binding affinity calculations indicated that 
hydrogen-bonding interactions played a key role 
between the TLR2 subfamily and malarial GPIs, simi- 
lar to lipopeptides. The identified key residues for 
binding of TLR2 subfamily members to malarial GPIs 
correlated with previous mutational and crystallogra- 
phy studies. The results obtained in the present study 
could provide new insights into understanding and reg- 
ulating the interaction between the TLR2 subfamily 
and malarial GPIs. 

Materials and methods 

Protein and ligand structure preparation 

Plasmodium falciparum Man4-GPI and the ,s7j-2 lyso GPI 
derivative were selected for molecular docking with TLR2- 
TLRl and TLR2-TLR6, respectively. The crystal coordi- 
nates of hTLR2-hTLRl-Pam3CSK4 (Protein Data Bank 
code: 2Z7X ) and mTLR2-mTLR6-Pam2CSK4 (Protein 
Data Bank code: 3A79 ) were used as the receptors after the 
ligands were removed using moe, version 2011.10 [45], The 
two-dimensional format of P. falciparum Man4-GPI and 
m-2 lyso GPI was drawn using chembiodraw ultra, ver- 
sion 12.0 (Cambridgesoft, Cambridge, MA, USA) and con- 
verted into the three-dimensional format by moe. The 
flexible alignment module in moe was used to align the 
structure of Man4-GPI to make it structurally similar to 
Pam3CSK4 and to make sn-2 lyso GPI structurally similar 
to Pam2CSK4, facilitating the accommodation of Man4- 
GPI and sn-2 lyso GPI into the same ligand-binding pock- 
ets of Pam3CSK4 and Pam2CSK4, respectively. 

Molecular docking 

The docking software gold, version 5.1 [46], from moe, ver- 
sion 2011.10, and Molegro virtual docker (mvd) software 
[47] were used, with Pam3CSK4 and Pam2CSK4 being used 
as reference ligands. gold uses a genetic algorithm for 
docking, and we selected the gold scoring function for the 
present study, mvd software uses an evolution algorithm 
and a piecewise linear potential scoring function. Before 
docking, the water molecules were removed to neutralize 
the receptor, ligands were protonated and partial charges 
were added by applying the MMFF94X forcefield [48] in 
MOE. The flexible small molecules Man4-GPI and sn-2 lyso 
GPI were docked into the binding sites of the macromole- 
cules hTLR2-hTLRl and mTLR2-mTLR6, respectively. 
Receptor coordinates of the crystal structure of hTLR2- 
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hTLRl-Pam3CSK4 and mTLR2-mTLR6-Pam2CSK4 com- 
plexes [17,18] were used to define the binding sites of 
Man4-GPI and sn-2 lyso GPI for docking simulations. The 
ligand binding sites were defined to detect the cavity within 
a sphere of 15-A radius around the coordinates of the 
ligands in the hTLR2-hTLRl-Pam3CSK4 and mTLR2- 
mTLR6-Pam2CSK4 complexes. The top ten docked 
conformations were saved among the 100 poses. Molecular 
interactions were observed via moe and poseview [49]. We 
used x-scoRE [37] to assess the binding affinity of the 
ligands. Images were generated using moe, chembiodraw, 
VMD [50] and chimera [51]. 

MD simulations 

The structures of hTLR2, hTLRl, mTLR2, mTLR6, 
hTLR2-hTLRl, mTLR2-mTLR6 and hTLR2-hTLRl- 
Pam3CSK4 in combination with docking studies were used 
to generate the hTLR2-hTLRl-Man4-GPI and mTLR2- 
mTLR6-.wi-2 lyso GPI complexes, which were subjected 
to 15 ns of MD simulations with yasara dynamics [52] 
using the amber03 force field [53]. Each protein was 
solvated in the centre of a cubic box with explicit [53] 
simple point charge water molecules, and periodic bound- 
ary conditions were applied in all directions. To replace 
the water molecules, appropriate Na^ or Cl~ counterions 
were used to neutralize each macromolecular system. The 
structures were energy minimized to remove bumps and 
fix the covalent geometry. In each step, energy minimiza- 
tion was performed by the steepest descent method in ya- 
sara dynamics. The simulations were performed at a 
constant temperature (300 °K) and pressure (1 bar) using 
the Berendsen coupling approach. For all simulations, a 
2-fs time step and 9-A nonbonded cut-off were used. The 
particle mesh Ewald method was used to treat long-range 
electrostatics, and bond lengths involving bonds to hydro- 
gen atoms were constrained using the SHAKE algorithm 
[54]. An additional 500-ps simulation and 15-ns produc- 
tion run were performed with snapshots collected every 
5 ps. The RMSD and RMSF were calculated using 
gromacs [55]. 

PCA 

PCA was used to examine the most substantial protein 
motions related to Ugand binding [56]. PCA allows the 
magnitude of a complicated protein system to be decreased, 
thus segregating the commanding modes of internal flux. 
The PCA of protein motions are computed as the eigenvec- 
tors of the mass-weighted covariance matrix of the back- 
bone atom's rearrangement. Protein coordinates of the 
TLRs were extracted from each time frame of the 5-ps MD 
trajectories for the protein complexes that were unbound 
and bound to GPIs. The g_covar gromacs tool was used 
to calculate and diagonalize the covariance matrix of back- 
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bone atoms to yield the eigenvalues and eigenvectors. The 
g_anaeig gromacs tool was used to analyze and plot the ei- 
genvectors to determine their dominant portion in protein 
motion. 
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